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Abstract 

' The scalar and vector self-energies obtained through QCD sum rules are introduced 

, in Quantum Hadrodynamics (QHD) equations. The results indicate that the ratios 

of coupling constants to respective meson mass have a very small dependence on 
ff^ ■ density. Then, an extended Brown-Rho scaling law is conjectured. 

o ■ 
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^ ' In recent years the density dependence of meson masses and meson-baryon 

coupling constants have received increasing interest [1-4]. Medium modifica- 
^ ■ tion of hadron properties can be seen as a manifestation of quark substructure. 

One of the principal reasons to study that topic is the success of Quantum 
Hadrodynamics (QHD) [5,6] in the relativistic description of nuclear phenom- 
ena. QHD has large success in the explanation of properties of infinite nuclear 
matter as well as finite nuclei. Several calculations of nuclear structure using 
QHD and its extensions were made obtaining a good agreement with experi- 
mental data [6,7]. 

QHD describes the N-N (nucleon-nucleon) interaction through mesons ex- 
change (tt, a, u, p, etc.). However, for the purposes of this work, only the 
simplest QHD model, namely QHD-I, which includes nucleons {ip) coupled 
with sigma [a) and omega (u) mesons, will be considered. In spite of the pion 
to be the principal component of the N-N interaction, it does not play any 
role in the model because nuclear matter is an isotropic system with parity 
conservation. Analyzing the model, it is possible to verify that the real part 
of the scalar term (cr-meson) is typically of the order of several hundred MeV 
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attractive while the real part of the time component of the vector term (u- 
meson) is typically of the order of several hundred MeV repulsive. However, 
the energies involved in problems of nuclear structure are only of a few tens 
of MeV. An energy of that order is obtained in QHD models due to a large 
cancellation between the scalar and vector pieces. In fact, scalar mesons pro- 
duce a long range attractive potential which saturate due to the presence of 
a repulsive short range potential produced by vector mesons. This process of 
saturation can be controlled through an appropriate choice of coupling con- 
stants. However, there are strong indications that such constants, as well as 
meson masses, change with the density [1,2]. So, the success of QHD, which 
uses fixed values for these amounts, must be explained. 

Recent results of the effective field theory (EFT) [6,8-12] have verified that 
QHD models are consistent with symmetries of Quantum Chromodynamics 
(QCD). Motivated by that, it will be very interesting to make a connection 
between QHD and QCD with the intention of studying coupling constants and 
meson masses in the medium. The natural way to make that connection take 
into account the nucleon self-energy. That idea was applied by Brockmann and 
Toki [2] to connect Relativistic Bruckner-Hartree-Fock Theory and Relativistic 
Hartree Approach. However, as it is known, the energies involved in nuclear 
matter problems are very lower than the energy scale of perturbative QCD. 
Thus, it is necessary to use a nonperturbative procedure. The indicated is 
the well-known method of QCD sum rules, that was introduced by Shifman, 
Vainshtein, and Zakharov in the late 1970's [13]. The QCD sum rule approach 
was successfully used in a lot of different problems, since hadron properties to 
particle decay. For the purposes of this work, the most interesting use of QCD 
sum rules is the possibility of describing properties of a nucleon in the medium. 
Indeed, QCD sum rules succeed to obtain the nucleon self-energy in terms of 
quarks and gluons degree of freedom [14-18]. On the other hand, self-energies 
obtained in QHD models are dependent on coupling constants and meson 
masses. The objective of this work is to study these hadronic parameters of a 
more fundamental point of view. 

To study meson-nucleon coupling constants and meson masses in nuclear mat- 
ter, it is first necessary to solve the nuclear ground state within Relativistic 
Hartree Approximation. The energy density in the QHD-I model can be writ- 
ten as [5] 

where M* is the effective nucleon mass in matter, kpis the Fermi momentum, 7 
is the spin-isospin degenerescence, gs and gy are, respectively, the cr-N and the 
a;-N coupling constants and and m„ are, respectively, the a- and a;-meson 
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masses. The self-energies E(s) and E(„) can be obtained through "tadpole" 
Feynman diagrams resulting in 

E(,) = M* - M = -4p. (2) 



and 



mi 



(3) 



Here M is the free nucleon mass with value M 939. MeV. The baryon and 
the scalar densities are, respectively, given by: 

PB = ^4 (4) 



and 



7M* 



kpEt - M*^ In 



kF + E*p 
M* 



(5) 



The self-energies (2) and (3) can also be obtained through QCD sum rules. The 
principal ingredient of this method is the time-ordered correlation function, 
defined by 

liapiq) =tjd^x e*«-^(0|T [vaix)r^m] |0) , (6) 



where |0) is the physical nonperturbative vacuum state and f]a{x) is an inter- 
polating field with the quantum numbers of a nucleon. In agreement with Ref. 
[18], we have for the proton field 



(7) 



where u and d are the up and down quark fields, a, b, c are color indices, 
T means transpose, C is the charge-conjugation matrix and t is an arbitrary 
real parameter. The OPE for the correlation function (6) can be generated 
through the following expansion for the quark propagator [19, 20] 



Sab{x) = {0\T[qa{x)qbm\0) 



. Sab /t 5ab > , 

^^^-T^(^?)-c + -- - , (8) 
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where {qq)vac is the free quark condensate, which can be determined from the 
Gell-Mann-Oakes-Renner relation, 

2m,{qq)vac = -mlf^ (l + 0(m^)) . (9) 



Here, m^r = 138 MeV is the pion mass, f^^ — 93 MeV is the pion decay constant 
and nig — {rriu + va^jl is the average quark mass. 

The same correlation function, Eq. (6), can be described through a phe- 
nomenological anzats inspired in the nucleon Green function. With the funda- 
mental state of the nuclear matter as vacuum and the interacting propagator 
in the phenomenological side, we obtain the following sum rules [18] 



^ U + 2t + + \E,{<iH),, ' ^^'^ 

^QCD _ 

where represent the Borel mass with value near ~ lGeV^[21], {q^q)pg — 
3pb/2 represents the quark density, {qq)pB is the in-medium quark condensate 



[22] and Eg = -S(,„) + j\^*2^ ^^^^ 191=270 MeV {i.e., approximately the 

Fermi momentum). Here, it was not take into account the continuous contri- 
bution. The expression (10), with t = —1, is a generalization of the lofFe's 
formula [21] to finite density. The exphcit dependence on the choice of the 
parameter t, that means, on the choice of the interpolating field, limits the 
significance of the loffc's formula to a qualitative role at best. A more realistic 
value for t can be obtained in order to do minimal the contributions com- 
ing from continuum model and high dimensional operators in OPE [24]. This 
better value will be considered later. Calculating the ratio of finite-density to 
zero-density sum rules, we obtain 



_ {qq)pn 
Eg5^^ 5 + 2t + 5t2 



(12) 

-PN . (13) 



Comparing expressions (2) and (12), we have for the scalar and vector sectors, 
respectively. 



9l 
ml 



{qq), 
{qq)i 



1 



(14) 
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and 



-M, 



N 



5 + 2t + 5t 



1 



(15) 



To analyze the expressions (12-15) is necessary to know the in-medium quark 
condensate, which can be written as 



{<l<l)pB = 1 - + • • ■ {m)vac , (16) 



where {qq)vac is given by Eq.(9) and the sigma term is estimated in Ref. [23] 
as being (Tn ~ 45 ± 10 MeV. The expression (16) is only been worth for low 
densities and considering those conditions we can write in a rude way the 
following approaches: ^ EqPB and pb ~ Ps- Therefore, the scalar sector 
is written as 



^QCD ctnMn 



g; gnMn Pb ctnM^ 



N 



TTT'lfn Ps rnlfl 
With the loffe's choice {t — —1), we obtain for the vector sector. 



(18) 



4 = ^- (20) 
ml mlf^ 

The principal conclusion obtained from equations (18) and (20) is that the 
ratios of coupling constants to respective meson mass are constants for low 
densities. Some values for these ratios are presented in table 1. The first hue 





all ml 


al/ml 


QHD-I 


3.029 10"^ 


2.222 10"^ 


QCDSR 


2.793 10-^ 


1.964 10-^ 



Table 1 

Ratios of coupling constants to respective meson mass. In the first line are presented 
the values used in Ref. [5] and in the second line the values obtained with equations 
(18) and (20) using lofFe's interpolating field [t = —1). 

of that table presents the ratios used in Ref. [5], which were adjusted so that 



5 



QHD-I reproduces the saturation properties of the nuclear matter. On the 
other hand, when the self-energies (17) and (19) and the ratios (18) and (20) 
(low density limit) are applied in Eq.(l), and we require the usual saturation 
condition for nuclear matter, namely e/ pb—M — 15.75 MeV at normal nuclear 
matter density {kp — 1-42 fm~^), we obtain the values presented in the second 
line of table 1. We can observe that the QCDSR results arc similar to the 
respective QHD-I values. However, this is not the most important result. The 
in-medium behavior of g1/m1 and gy/ml, named scalar and vector ratios, 
respectively, are the principal subject to be studied here. 

Despite the fact of we have obtained the equilibrium properties of nuclear 
matter, the effects despised in the low density approximation must be included. 
For that, it is necessary to abandon Eq.(16) for the condensate and to use 
the familiar Nambu-Jona-Lasinio (NJL) model, which give us the following 
expression [22] for the quark condensate: 



The Fermi momentum of quarks, kpf^^^ is related to the Fermi momentum of 
nucleons through the density: Pg — 3pB, where Pg is the quark density. Here 
iVc is the color number, AjyjL is a cutoff, Gnjl is a coupling constant, nig 
is the current quark mass and Mg is the constituent quark mass, which is 
dynamically generated by a partial restoration of chiral symmetry. With an 
appropriate choice of these parameters the NJL model gives us good results 
for the TT-meson mass, f^,, constituent quark mass, condensate, etc. A good 
choice is A^vjl = 900 McV and G^jl = 3.54/A^;^, for which the saturation 
point of nuclear matter is reproduced. The in-medium behavior of the relevant 
quantities g1/m1 and gy/ml, are presented on figure 1. The scalar [Eq. (14)] 
and the vector [Eq.(15)] ratios are represented by the solid and the dashed 
lines, respectively. In these calculation, a more realistic nucleon interpolating 
field (i = -1.1) [24] has been used. 

In summary, the expressions (17-20) are valid just for the low density region. 
However, some insight can be obtained from these first results: the ratios 
of coupling constants to respective meson masses are constants in the low 
density limit. Actually, the principal conclusion obtained in this work can be 
summarized by figure 1, where the scalar and the vector ratios are presented 
as functions of the density. We can observe that they have a small dependence 




(21) 



and 



Mg^ruq- 2GNjL{qq) 



(22) 
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Fig. 1. The ratios g1/m1 (solid line) and g^jvn^ (dashed line) calculated with 
Eqs.(14) and (15), respectively 

on density and the vector ratio decrease when the density is enlarged. These 
results corroborate with the observation made in Ref.[4], where the nucleon 
flow probing higher density requires that gv/ray be independent of density at 
low densities and decrease slightly at high densities. The results here obtained 
indicate that {gs/'ms)vac ~ {.gs/i^s)pB "^^ can write 



{9s 
{9s)x 



'PB 



(23) 



The same result is valid for the vector ratio. In other words, couphng constants 
must change with the same rate of the respective meson mass, in order to keep 
the ratios unaffected. Therefore, the fact of the ratios remain constants can 
give us an explanation of how QHD models, that use parameters independent 
of density, have success to explain the bulk properties of nuclear matter. The 
same is not true for finite nuclei, where the scalar and vector ratios do not 
appear explicitly on equations. Thus, the dependence on density of coupling 
constants and meson masses must be important in these cases. 

Equation (23) and its equivalent for the vector sector, permit one to extend the 
Brown-Rho (BR) scaling law [1] to include coupling constants in the following 
way: 

Mk^ll^<^<^9l^9:^ .24) 
Mn U ms my gs gv' 
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where the asterisks denote in-medium quantities. 

Finally, it is necessary to say that these results are valid at low densities. For 
nuclear matter in extreme conditions, the ratios should change [4] and QHD- 
I must present bad results. Besides, in these calculations the contributions 
of gluon condensates and other high terms were not included in the OPE. 
However, we know that the contributions of these terms for the sum rule are 
very small. The use of an extended version of the sum rules (10) and (11) as 
well as the importance of scaling for finite nuclei, are left for future works. 
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